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ABSTRACT 

Context. The mass-loss rate is a key parameter of hot, massive stars. Small-scale inhomogeneities (clumping) in the winds of these 
stars are conventionally included in spectral analyses by assuming optically thin clumps, a void inter-clump medium, and a smooth 
velocity field. To reconcile investigations of different diagnostics (in particular, unsaturated UV resonance lines vs. H„/radio emis- 
sion) within such models, a highly clumped wind with very low mass-loss rates needs to be invoked, where the resonance lines seem 
to indicate rates an order of magnitude (or even more) lower than previously accepted values. If found to be realistic, this would 
challenge the radiative line-driven wind theory and have dramatic consequences for the evolution of massive stars. 
Aims. We investigate basic properties of the formation of resonance lines in small-scale inhomogeneous hot star winds with non- 
monotonic velocity fields. 

Methods. We study inhomogeneous wind structures by means of 2D stochastic and pseudo-2D radiation-hydrodynamic wind mod- 
els, constructed by assembling ID snapshots in radially independent slices. A Monte-Carlo radiative transfer code, which treats the 
resonance line formation in an axially symmetric spherical wind (without resorting to the Sobolev approximation), is presented and 
used to produce synthetic line spectra. 

Results. The optically thin clumping limit is only valid for very weak lines. The detailed density structure, the inter-clump 
medium, and the non-monotonic velocity field are all important for the line formation. We confirm previous findings that radiation- 
hydrodynamic wind models reproduce observed characteristics of strong lines (e.g., the black troughs) without applying the highly 
supersonic 'microturbulence' needed in smooth models. For intermediate strong lines, the velocity spans of the clumps are of central 
importance. Current radiation-hydrodynamic models predict spans that are too large to reproduce observed profiles unless a very low 
mass-loss rate is invoked. By simulating lower spans in 2D stochastic models, the profile strengths become drastically reduced, and 
are consistent with higher mass-loss rates. To simultaneously meet the constraints from strong lines, the inter-clump medium must be 
non-void. A first comparison to the observed Phosphorus V doublet in the 06 supergiant A Cep confirms that line profiles calculated 
from a stochastic 2D model reproduce observations with a mass-loss rate approximately ten times higher than that derived from the 
same lines but assuming optically thin clumping. Tentatively this may resolve discrepancies between theoretical predictions, evolu- 
tionary constraints, and recent derived mass-loss rates, and suggests a re-investigation of the clump structure predicted by current 
radiation-hydrodynamic models. 

Key words, stars: early-type - stars: mass-loss - radiative transfer - line: formation - hydrodynamics - instabilities 



O ■ 1 ■ Introduction Il995t iDessart & Owockill2005l) . where the line-driven (or line- 

L" ' deshadowing) instability causes a small-scale, inhomogeneous 

Mass loss through supersonic stellar winds is pivotal for the wind fa both densky and velodty 



physical understanding of hot, massive stars and their surround 
ings. A change of only a factor of two in t he mass-loss rate ha s Direct observational evidence of a small-scale, clumped stel 



x 

■ a dramatic effect on massive star evolution (M evnet et al.lll994l) . lar wind has, for O-star s so far only been given for two objects. 

Win ds from these stars are described by the lin e-driven wind the- I Pup and HD 93 129 A ( lEversberg et al.|[T998t iLepine &~M offat 

ory (ICastor et al.ll 19751: iPauldrach et al.ll 1986k . which tradition- 2008). Much indirect evidence, however, has arisen from quan- 

ally assumes the wind to be stationary, spherically symmetric, titative spectroscopy, where the standard way of deriving mass- 

and homogeneou s. Despite this theory's apparent success (e.g., loss rates from observations nowadays is via line-blanketed, 

IVink et"dTl2000b . evidence for an inhomogeneous and time- non-LTE (LTE: local thermodynamic equilibrium) model 



at- 

dependent wind has over the past years accumulated, recently mospheres that include a treatment of both the photosphere 

summarized in the proceedings from t he workshop 'Clumping and the wind. W ind clumping has been include d in such codes 

in hot star winds' dHamann et al.ll2008l) and in a gener al review (e.g., CMFGEN dHillier & Millerlll998h . PoWR dGrafener et al.l 

of mass loss from hot, massive stars dPuls et"aT]|2008bl) . 2002), FASTWIND dPuls et alJl2005h ') bv assuming statistically 

That line-driven wi nds should be intr i nsica lly unstable was distributed optically thin density clumps and a void inter-clump 

already pointed out bv lLucv & Solomonl d 19701) . and was later medium, while keeping the smooth velocity law. The major 

confirmed first by linear stability analyses and then by direct, result from this methodology is that any mass-loss rate de- 

radiat ion-hydrodynamic model i ng of the time-dependent wind rived from smooth models and density-squared diagnostics (H a , 

(e.g., lOwocki & Rvbickil [l984t lOwocki et al.l IT988: Feldmeier infra-red and radio emission) needs to be scaled down by the 
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square root of the clumping factor (which describes the over 
density of the clumps as compared to t h e mea n density, see 
Sect. 2~2l l. For example. ICrowther et all d2002l) . iBouret et al.l 
(2003), and Bo uret etall (120051) have concluded that a reduction 
of 'smooth' mass-loss rates by factors 3 ... 7 might be neces- 
sary. Furthermore, from a combined optical /IR/radio analys is of 
a sample of Galactic O-giants/supergiants. lPuls et alJ (120061) de- 
rived upper limits on observed rates that were factors of 2 ... 3 
lower than previous H a estimates based on a smooth wind. 

On the other hand, the strength of UV resonance lines ('P 
Cygni lines') in hot star winds depends linearly on the den- 
sity and is therefore not believed to be directly affected by op- 
tically thin clumping. By using the Sob olev with exact inte- 
gration technique (SEI; cf. lLamers et alj 1987) o n the unsat- 
urated Phosphorus V (PV) lines. iFullerton et all d2006l) for a 
large number of Galactic O-stars derived rates that were factors 
of 10 ... 100 lower than corresponding smooth H ff /radio values 
(provided PV is the dominant ion in spectral classes 04 to 07). 
Such large revisions would conflict with the radiative line-driven 
wind theory and have dramatic consequences for the evolution 
of, and the feedback from, massive stars (cf. ISmith & Owockil 
2006: lHirschill2008h . Indeed, a puzzling picture has emerged, 
and it appears necessary to ask whether the present treatment of 
wind clumping is sufficient. Particularly the assumptions of op- 
tically thin clumps, a void inter-clump medium, and a smooth 
velocity field may not be adequate to infer proper rates under 
certain conditions. 



Optically thin vs. optically t hick clumps. lOskinoya et all d2007l) 
used a porosity formalism dFeldmeier et al. 12003 : Owocki et alJ 
|2004|) to scale the opacity from smooth models and investigate 
impacts from optically thick clumps on the line profiles of £ Pup. 
Due to a reduction in the effective opacity, the authors were able 
to reproduce the PV lines without relying on a (very) low mass- 
loss rate, while simultaneously fitting the o ptically thin H ? line. 
This formalism, however, was criticized bv lOwockl (12008) who 
argued that the original porosity concept had been developed for 
continuum processes, and that line transitions rather should de- 
pend on the non-monotonic velocity field seen in hydrodynamic 
simulations. Proposing a simplified analytic description to ac- 
count for this velocity -porosity, or 'vorosity', he showed how 
also this effect may reduce the effective opacity. 

In this first paper we attempt to clarify the most important 
concepts by conducting a detailed investigation on the synthe- 
sis of UV resonance lines from inhomogeneous two-dimensional 
(2D) winds. We create both pseudo-2D, radiation-hydrodynamic 
wind models and 2D, stochastic wind models, and produce 
synthetic line profiles via Monte-Carlo radiative transfer cal- 
culations. We account for and analyze the effects from a wind 
clumped in both density and velocity as well as the effects from 
a non-void inter-clump medium. Especially we focus on lines 
with intermediate line strengths, comparing the behavior of these 
lines with the behavior of both optically thin lines and saturated 
lines. Follow-up studies will include a treatment of emission 
lines (e.g., H„) and an extension to 3D, and the development of 
simplified approaches to incorporate effects into non-LTE mod- 
els. 

In Sect. [2] we describe the wind models and in Sect. [3] the 
Monte-Carlo radiative transfer code. First results from 2D inho- 
mogeneous winds are presented in Sect. [4] and an extensive pa- 
rameter study is carried out in Sect. [5] We discuss some aspects 
of the interpretations of these results and perform a first compar- 



ison to observations in Sect. [6] and summarize our findings and 
outline future work in Sect. [7] 



2. Wind models 

For wind models, we use customary spherical coordinates 
(r, 0, <5) with r the radial coordinate, © the polar angle, and <J the 
azimuthal angle. We assume spherical symmetry in ID models 
and symmetry in <t> in 2D models. In all 2D models is sliced 
into Nq equally sized slices, giving a lateral scale of coherence 
(or an opening angle) 180/Af© degrees. This 2D approximation 
is discussed in Sect. 16.41 Below we describe the model types pri- 
marily used in the present analysis; two are of stochastic nature 
and two are of radiation-hydrodynamic nature. 



2.1. Radiation-hydrodynamic wind models 

We use the time-depen dent, radia t ion-hy drodynamic (hereafter 
RH) wind models from lPuls etal] d!993l hereafter 'POF' ), cal- 
culated by S. Owocki, and from Feldmeier et al] dl997l here- 
after 'FPP'), and the reader is referred to these papers for de- 
tails. Here we summarize a few important aspects. POF as- 
sume a ID, spherically symmetric outflow, and circumvent a 
detailed treatment of the wind energy equation by assuming 
an isothermal flow. Perturbations are triggered by photospheric 
sound waves. The wind consists of 800 radial points, extend- 
ing to roughly 5 stellar radii. FPP also assume a ID, spherically 
symmetric outflow, but include a treatment of the energy equa- 
tion. Perturbations are triggered either by photospheric sound 
waves or by Langevin perturbations that mimic photospheric tur- 
bulence. The wind consists of 4000 radial points, extending to 
roughly 30 stellar radii. Tests have shown that the FPP winds 
yield similar results for both flavors of perturbations, and, for 
simplicity, we therefore use only the results of the turbulence 
model. 

Due to the computational cost of obtaining the line force, 
only initial attempts to 2D RH si mulations have been carried 
out (Dessart & Owocki 200l l200l . These authors first used a 
strictly radial line force, yielding a complete lateral incoherent 
structure due to Rayleigh-Taylor or thin-shell instabilities, and in 
the follow-up study uses a restricted 3-ray approach to approxi- 
mate the lateral line drag, yielding a larger lateral coherence but 
lacking quantitative results. Therefore, and because of the gen- 
eral dominance of the radial component in the radiative driving, 
we create fragmented 2D wind models from our ID RH ones 
by assembling snapshots in the direction, assuming indepen- 
dence between each slice consisting of a pure radial flow. After 
the polar angle has been sliced into N@ equally sized slices, one 
random snapshot is selected to represent each slice. This method 
for creating more-D models fr om ID ones is essent i ally th e same 
as the 'patch method' used by iDessart & Owockil d2002l) . when 
synthesizing e mission lines fo r Wolf- Rayet stars, and the method 
used by, e.g., lOskinova et al.l d2004l) . when synthesizing X-ray 
line emission from stochastic wind models. Fig.[T]displays typi- 
cal velocity and density structures from this type of 2D model. 



2.2. Stochastic wind models 

We also study clumpy wind structures created by means of dis- 
torting a smooth, stationary, and spherically symmetric wind via 
stochastic procedures. This allows us to investigate the impacts 
from, and to set constraints on, different key parameters with- 
out being limited by the values predicted by the RH Simula- 
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Fig. 1. Left panel: Density contour plots of one stochastic (upper plot) and one RH (FPR lower plot) model. The Cartesian coordinate 
Z is on the abscissa and X is on the ordinate. Right panel: Density and velocity structures of one slice in one stochastic (upper) and 
one RH (FPR lower) model. Over densities are marked with filled dots. For model parameters and details, see Sect. 12.21 



tions. For the underlying smooth winds we adopt a standard /3 
velocity law vp(r) = (1 - b/rf. Here and throughout the paper, 
we measure all velocities in units of the terminal velocity, v M , 
and all distances and length scales in units of the stellar radius, 
R+. b is given by v(r — 1) = v m j n , the velocity at the base of 
the wind. v m ; n = 0.01 is assumed, roughly corresponding to the 
sound speed. For a given M, the homogeneous density structure 
then follows directly from the equation of continuity. We choose 
f3 = 1, which is appropriate for a standard O-star wind and al- 
lows us to derive simple analytic expressions for wind masses 
and flight times. 



A model clumped in density. First we consider a two compo- 
nent density structure consisting of clumps and a rarefied inter- 
clump medium (hereafter ICM), but keep the ft = 1 velocity law. 
Clumps are released randomly in radial direction at the inner 
boundary, independently from each slice. The release in radial 
direction means that a given clump stays within the same slice 
during its propagation through the wind. The average time inter- 
val between the release of two clumps is 6t, which here and in 
the following is expressed in units of the wind's dynamic time 
scale f dyn = ^*/voo. 



The average distance between clumps thus is VpSt, i.e. 
clumps are spatially closer in the inner wind than in the outer 
wind, and for example 5t = 0.5 (in fd yn ) gives an average clump 
separation of 0.5 (in 7?*) at the point where v = 1 (in Voo). We 
further assume that the clumps preserve mass and lateral angle 
when propagating outwards, and that the underlying model's to- 
tal wind mass is conserved within every s lice. This radial clump 
distribution is the same as the one used bv lOskinova et al] 12006) 
when simulati ng X-ray emission from O-stars, but differs from 
the one used bv lOskinova et al](l2007l) when investigating poros- 



ity effects on resonance lines (see discussion in Sect. 16.5b . The 
radial clump widths are here calculated from the actual wind 
geometry and clump distribution by assuming a volume filling 



factor f y , defined as the fractional volume of the dense gafl A 
related quantity is the clumping factor 



/c1 = <p> 2 ' 



(i) 



as defined bv lOwocki etal] d!988l) . where angle brackets denote 
temporal averages. Identifying temporal with spatia l averages 
one m ay write for a two component medium (cf. lAbbott et al.l 
[19811) 



/v+(l-/v)*? 
[/ V + (1-/ V )X 1C ] 2 



with 



Pel 



(2) 



(3) 



the ratio of low- to high-density gas (subscript ic denotes inter- 
clump and cl denotes clump). For a void (jq c = 0) ICM, p c i/(p) 
= f~ l - / c i, i.e, / c i then describes the over density of the clumps 
as compared to the mean density. 

A model clumped in density and velocity. Next we consider 
also a non-monotonic velocity law, using the spatial distribution 
and widths of the clumps described in the previous paragraph. 
The RH simulations indicate that, generally, strong shocks sepa- 
rate denser and slower material from rarefied regions with higher 
velocities. Building on this basic result, we now modify the ve- 
locity fields in our stochastic models by adding a random per- 
turbation to the local Vp value prior to the starting point of each 
clump, so that the new velocity becomes v pre . A 'jump velocity' 
is thereafter determined by a random subtraction from vp, now 
using the added perturbation as the maximum subtraction. That 
is, 



"pre 



vp + vj x 2Ri 



v post 



vp - Vj x 2RiR 2 , 



(4) 



1 We here notice that / v is normalized to the total volume, i.e., f y = 
... 1. In some literature /„ is identified with the straight volume ratio 
Vd/Vic which then implicitly assumes that V c i <k Vj c . 
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Table 1. Basic parameters defining a stochastic wind model 
clumped in density and with a non-monotonic velocity field. 



Name 



Parameter Considered range 



Volume filling factor f v 

Average time interval be- St 
tween release of clumps 
ICM density parameter, Eq.|3] x K 

Velocity span of clump Sv 

Parameter determining the vj 
jump velocity 



/ v =0.01... 1.0 
6t[t iy n]= 0.05... 1.5 

x ic =0...0.1 
Sv/Svp = -10.0... 1.0 
v-Jvp =0.01... 0.15 




that treats resonance line formation in a spherical and axially 
symmetric wind using an 'exact' formulation (e.g., without re- 
sorting to the Sobolev approximation). The restriction to 2D is 
of course a shortage, but has certain geometrical and computa- 
tional advantages and should be sufficient for the study of gen- 
eral properties, as discussed in Sect. 16.41 A thorough description 
and verification of the code can be found in AppendixlAl 

Photons are released from the lower boundary (the photo- 
sphere) and each path is followed until the photon has either 
left the wind or been backscattered into the photosphere. Basic 
assumptions are a line-free continuum with no limb darkening 
emitted at the lower boundary, no continuum absorption in the 
wind, pure scattering lines, instantaneous re-emission, and no 
overlapping lines (i.e., singlets). These simplifying assumptions, 
except for doublet formation, are all believed to be of minor 
importance to the basic problem. By the restriction to singlet 
line formation we avoid confusion between effects on the line 
profiles caused by line overlaps and by other important param- 
eters, but on the other hand it also prevents a direct comparison 
to observations for many cases (but see Sect. 16.6b . A consistent 
treatment of doublet formation will be included in the follow-up 
study. 



Fig. 2. Non-monotonic velocity field and corresponding param- 
eters in a stochastic model. 



where R\ and R2 are two random numbers in the interval to 1. 
p 0s t is the jump velocity as determined by the parameter 



v, 



pre 

Vj. By multiplying R\ by two, we make sure that the mean per- 
turbation at the 'pre' point is vj, and R2 allows for an asymmetry 
about vp (see Fig. |2j. The clump is assumed to start at v post , and 
its velocity span is set by assuming a value for dv/dvp, where 
Sv is the velocity span of the clump and 5vp the corresponding 
quantity for the same clump with a smooth velocity law (see 
Fig. |2}. Inspection of our RH models suggests that velocity gra- 
dients within density enhancements primarily are negative (see 
also Sect. 16.3b . and negative gradients are also adopted in most 
of our stochastic models. Finally we assume a constant velocity 
gradient through the ICM. 

Overall, the above treatment provides a phenomenological 
description of the non-monotonic velocity field seen in RH 
simulations. T he description differs from the one suggested by 
lOwockil ([2008), who uses only one parameter to characterize the 
velocity field (whereas we have two). Our new formulation is 
motivated by both observational and modeling constraints from 
strong and intermediate lines, as discussed in Sect. 16.51 

The basic parameters defining a stochastic model are listed 
in Table Q] Fig. Q] (right panel) shows the density and velocity 
structures of one slice in a stochastic model, with density pa- 
rameters / v = 0.1, St = 1.0, Xi c = 0.005, and velocity parameters 
vj = Q.\5vp and 6v - —dvp. Clump positions have been high- 
lighted with filled dots and a comparison to a RH model (FPP) 
is given. In the RH model, we have identified clump positions 
by highlighting all density points with values higher than the 
corresponding smooth model. The left panel shows the density 
contours of the same models, where, for clarity, only the wind to 
r - 5 is displayed. 



3. Radiative transfer 

To compute synthetic line profiles from the wind models, we 
have developed a Monte-Carlo radiative transfer code (MC-2D) 



4. First results from 2D inhomogeneous winds 

Throughout this section we assume a thermal velocity, v t = 
0.005 (in units of v«, and ~ 10 kms -1 , appropriate for a stan- 
dard O-star wind), and apply no micro turbulence. After a brief 
discussion on the impact of the observer's position and opening 
angles, we concentrate on investigating the formation of strong, 
intermediate, and weak lines. In our definition, an intermediate 
line is characterized by a line strength^ Kq =5.0 chosen such as 
to almost precisely reach the saturation limit in a smooth model 
(cf.Fig.E). 

By investigating these different line types, we account for the 
tight constraints that exist for each flavor: i) weak lines should 
be independent of density-clumping properties as long as the 
clumps remain optically thin, ii) for intermediate lines either 
smooth models overestimate the profile strengths or mass-loss 
rates are lower than previously thought (e.g. the PV problem, see 
Sect. [TJ, and iii) strong saturated lines are clearly present in hot 
star UV spectra, and observed features need to be reproduced, 
such as high velocity (> v^) absorption, the black absorption 
trough, and the reduction of re-emitted flux blueward of the line 
center. 



4.1. Observer's position and opening angles 

The observed spectrum as calculated from a 2D wind structure 
depends on the observer's placement relative to the star (see 
Appendix lAli. As it turns out, however, this dependence is rel- 
atively weak in both the stochastic and the RH models (the latter 
is demonstrated in the upper panel of Fig. [3}. Tests have shown 
that the variability of the line profile's emission part is insignif- 
icant. The variability of the absorption part may be detectable, 
at least near the blue edge, but is still insignificant for the in- 
tegrated profile strength; the equivalent width of the absorption 
part is almost independent of the observer's position. Also the 
opening angle, 180° /Ms, primarily has a smoothing effect on the 
profiles. In Fig. [3] prominent discrete absorption features appear 



2 with kq proportional to the product of mass-loss rate and abundance 
of the considered ion, see AppendixlAl 
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Fig. 3. Synthetic line profiles calculated from 2D RH models. 
The abscissa is the dimensionless frequency x (Eq. IA.1 11 1, nor- 
malized to the terminal velocity, and the ordinate is the flux 
normalized to the continuum. Upper panel: Profiles from POF 
models with kq = 5.0. The upper plot displays profiles for an 
observer placed at the o b. s angles as labeled in the figure and a 
profile averaged over all N@ = 30 angles. The lower plot displays 
averaged profiles for three different Nq. Lower panel: Averaged 
profiles from FPP and POF models with N@ = 30, and with 
kq = 100.0 (upper) kq = 5.0 (middle), and kq = 0.05 (lower). For 
comparison, ID, smooth profiles with the same values of kq are 
shown as well. 



near the blue edge in the model with N@ — 1 (spherical sym- 
metry), but are smoothed out in the 'broken-shell' models with 
Nq = 30 and 60. The equivalent widths of the absorption parts 
are approximately equal for all three models. 

Because our main interest here is the general behavior of 
the line profiles, we choose to work only with Nq = 30 and 
profiles averaged over all observer angles from here on. Working 
with averaged line profiles has great computational advantages, 
because roughly a factor of N@ fewer photons are needed. 

4.2. Radiation-hydrodynamic models 

Fig. [3] (lower panel) shows line profiles from FPP and POF hy- 
drodynamical models. For the strong lines, the constraints stated 



in the beginning of this section are reproduced without adopting 
a highly supersonic and artificial microturbulence. These fea- 
tures arise because of the multiple resonance zones in a non- 
monotonic velocity field, and are present in spherically symmet- 
ric RH profiles as well (see POF for a comprehensive discus- 
sion); the main difference between ID and 2D is a smoothing 
effect, partly stemming from averaging over all observer angles 
(see above). The absorption at velocities higher than the terminal 
is stronger in FPP than in POF, due to both a higher velocity dis- 
persion and a larger extent of the wind (r max ~ 30 as compared 
to r max ~ 5, see Sect. 12. Ik more overdense regions are encoun- 
tered in the outermost wind, which (because of the flatness of 
the velocity field) leads to an increased probability to absorb at 
almost the same velocities. 

For the intermediate lines, we again see the qualitative fea- 
tures of the strong lines, though less prominent. As compared 
to smooth models, a minor absorption reduction is present at 
velocities lower than the terminal, but compensated by the blue 
edge smoothing. Therefore the equivalent width of the line pro- 
file's absorption part in the FPP model is approximately equal 
to that of the smooth model, whereas in the POF model it is re- 
du ced by ~ 10%. This minor reduction agrees with that found 
bv lOwockil (12008). and is not strong enough to explain the ob- 
servations without having to invoke a very low mass-loss rate. 

For the weak lines, the absorption part is marginally stronger 
than from a smooth, ID model. 

4.3. Stochastic models 

In this subsection we use a 'default' 2D, stochastic model with 
parameters as specified in Table [2] By comparing this model 
to models in which one or more parameters are changed, we 
demonstrate key effects in the behavior of the line profiles. 

Strong lines. For strong lines, the line profiles from the default 
model reproduce the observational constraints described in the 
first paragraph of this section. As in the RH models, we apply 
no microturbulence. Fig.|4](left panels) demonstrates the impor- 
tance of the ICM in the default model; the absorption part of a 
very strong line is not saturated when x lc = 0. That is, with a 
void ICM we will, regardless of the opacity, always have line 
photons escaping their resonance zones without ever interacting 
with any matter, thereby de-saturating the line. This ICM find- 
ing agrees with that of Zsar go et alj (12008), who point out that 
a non-void ICM is crucial for the formation of highly ionized 
species such as O VI. We also notice that 5v - —Svp (used in the 
default model) does not permit clumps to have velocities higher 
than the local Vp value, preventing absorption at velocities higher 
than the terminal one when the ICM is void. 

Intermediate lines. For intermediate lines, the line profiles 
from the default model display the main observational require- 
ment if to avoid a drastic reduction in 'smooth' mass-loss ratefl 
namely a strong absorption reduction as compared to a smooth 
model. The left panels of Fig. [4] show how the integrated profile 
strength of the default model with Kq =5.0 roughly corresponds 
to that of a smooth model having kq = 0.5, i.e., the smooth model 
would result in a mass-loss rate (as estimated from the integrated 
profile strength) ten times lower than the clumped model. The 
figure also illustrates how the main effect is on the absorption 

3 Recall that / v = 0.25 — > f a ~ 4, which implies M = M sm0 oth/2, if 
f c \ were derived from p 2 -diagnostics assuming optically thin clumps. 
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Table 2. Primary stochastic wind models and parameters 



Model name 


/v 


St [t dyn ] 


Xic 








r b 

' ext 


Default 


0.25 


0.5 


0.0025 


-1.0 


0.15 


1.3 


-25 


RHcopy 


0.1 


0.5 


0.005 


-10.0 


0.15 


1.3 


~5 


Obsl 


0.11 


0.5,4.0 c 


0.005,0.0025 c 


-1.0 


0.15 


1.02 


-25 



a Radial onset of clumping. b Radial extent of wind. 

c Left value inside the radius corresponding to vp = 0.6, right value outside. 
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Fig. 4. Left panels: Solid lines 
display total line profiles and 
the absorption part for the 
default stochastic model (see 
Table 0, with k = 1000 
(upper), kq = 5.0 (mid- 
dle), and k = 0.05 (lower). 
Dotted lines display smooth 
models with k = 5.0 and 
kq = 0.5 (middle), and kq = 
0.05 (lower). Dashed/dashed- 
dotted lines with modifica- 
tions from the default model 
as labeled in the figure. 
Right panels: Same as the 
left panels, but for POF 
(dashed lines) and RHcopy 
(solid lines). Dashed-dotted 
lines with modifications from 
RHcopy as labeled in the fig- 
ure. 



part of the line profile. In addition to the reduction in profile 
strength, the profile shapes of the absorption parts are noticeably 
different for the default and smooth models (the shapes of the re- 
emission parts, not shown here, are similar for the two models). 
We further discuss the shapes of the profiles in Sect. 16.11 The 
dramatic reduction in integrated profile strength occurs because 
of large velocity gaps between the clumps, in which the wind is 
unable to absorb (at this opacity the ICM may not 'fill in' these 
gaps with absorbing material). 

We have identified \Sv\ as a critical parameter for the forma- 
tion of intermediate lines. The importance of the velocity spans 
of the clumps is well illustrated by the absorption part profiles 
in Fig. |4] (lower-left panel, middle plot). The absorption is much 
stronger in the comparison model with 8v = -55vp than in the 
default model with 5v = — 5vp, because the former model covers 
more of the total velocity space within the clumps, thereby clos- 
ing the gaps between the clumps. Consequently the wind may, 
on average, absorb at many more wavelengths. 

In principle, however, this effect is counteracted by a de- 
crease in the clump's optical depths, because of the now higher 
velocity gradients (\dv/8\>p\ > 1). Consider the radial Sobolev 
optical depth (proportional to p/\dv/dr\, see Appendix [Al in a 
stochastic wind model. As compared to a smooth model, the 
density inside a clump is enhanced by a factor of f~ l (assum- 



ing a negligible ICM), but also the velocity gradient is enhanced 
by a factor of \5v/6vp\. Thus we may write for the radial Sobolev 
optical depth inside a clump, 

^"Sob.sm ^0 , r -. 

TSob ~ 1 ~ , (5) 

fy\Sv/dVfi\ Vfify\Sv/SVfi\ 

where 'sm' indicates a quantity from a smooth wind, and the 
expression to the right is valid for an underlying [3-1 veloc- 
ity law. From Eq. [5] we see how the effects on the optical depth 
from the increased density (/ v = 0.25) and the increased velocity 
gradients (\dv/Svp\ = 5) almost cancel each other in this exam- 
ple. Thus, the clumps are still optically thick for the intermediate 
line (ko = 5), which means that the larger coverage of the total 
velocity space 'wins', and the net effect becomes an increase in 
absorption (as seen in Fig. |4] lower-left panel, middle plot). This 
will be true as long as not f v \Sv/6vp\ » 1, which is never the 
case in the parameter range considered here. 

Finally, the prominent absorption dip toward the blue edge 
in the default model turns out to be a quite general feature of our 
stochastic models, and is discussed in Sects. 15. 1 l and !6.2l 

Weak lines. The statistical treatment of density clumping in- 
cluded in atmospheric codes such as CMFGEN, PoWR, and 
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FASTWIND is valid for optically thin clumps and a negligible 
ICM, and gives no direct effect on resonance lines scaling lin- 
early with density. Here we test this prediction using detailed ra- 
diative transfeiQ. Our default model recovers the smooth results 
when kq = 0.05 (Fig. [4] left panels), confirming the expected 
behavior. However, from calculating spectra using different val- 
ues of kq, we have found that significant deviations from smooth 
models occur for the default model already before a-o reaches 
unity. This occurs because the clumps start to become optically 
thick, which may again be understood by considering the radial 
Sobolev optical depth (Eq.|5]). With / v < 0.25 and k > 0.25, one 
finds Tsob ^ 1 -0. 

4.4. Comparison between stochastic and 
radiation-hydrodynamic models 

Our stochastic wind models have been constructed to contain all 
essential ingredients of the RH models. Therefore they should 
also reproduce the RH results, at least qualitatively, if a suitable 
parameter set is chosen. To test this we used the POF model. In 
this model, the clumping factor increases drastically at r ~ 1.3, 
from / c i ~ 1 .0 to / c i ~ 10, after which it stays basically constant. 
The average clump separation in the outer wind is roughly half a 
stellar radius. Important for the velocity field is that the velocity 
spans of the clumps are generally larger than corresponding 
spans', i.e., \8v/Svp\ > 1 (this is the case in FPP as well), a char- 
acteristic behavior that primarily affects the intermediate lines 
(details will be discussed in Sect. I6.31 l. Finally, a suitable Vj can 
be assigned from the position of the blue edge in a strong line 
calculated from POF. Table [2] (entry RHcopy) summarizes all 
parameters used to create this stochastic, 'pseudo-RH' model. 
Fig.|5]displays one slice of the velocity and density structures in 
the POF and RHcopy models, and Fig.|4](right panels) displays 
the line profiles. 

The line profiles of POF are matched reasonably well by 
RHcopy. The intermediate lines again demonstrate the impor- 
tance of the velocity spans of the clumps; for an alternative 



Fig. 6. Left: Schematic of Av, the velocity gap between two sub- 
sequent clumps, made up by velocities not covered by any of the 
clumps. 6v is the velocity span of a clump and v t the thermal ve- 
locity. Right: The effective escape ratio rj (Eq. [7} as a function 
of /3 = 1 velocity, for the parameters of the default model (see 
Table©. 

model with 6v = -8vp, there is much less absorption in the 
stochastic model than in POF, i.e., we encounter the same effect 
as discussed in the previous subsection. We conclude that in RH 
models it is the large velocity spans inside the density enhance- 
ments that prevent a reduction in profile strength (as compared 
to smooth models) for intermediate lines. 

5. Parameter study 

Having established basic properties, we now use our stochastic 
models to analyze the influence from different key parameters 
in more detail. First, however, we introduce a quantity that turns 
out to be particularly useful for our later discussion. 

5.1. The effective escape ratio 

For the important intermediate lines, it is reasonable to assume 
that the clumps are optically thick and the ICM negligible (see 
Sect. I4.3I and the next paragraph). Under these assumptions, a 
decisive quantity for photon absorption will be the velocity gap 
not covered by the clumps, as compared to the thermal velocity 
(the latter determining the width of the resonance zone in which 
the photon may interact with the wind material). This is illus- 
trated in the left panel of Fig. [6] and we shall call this quantity 
the 'effective escape ratio' 



Av 



n = 



(6) 



4 The indirect effect through the feedback on the occupation numbers 
is not included, because in this section we assume constant ionization. 



where Av is the velocity gap between two subsequent clumps, 
made up by all velocities not covered by any of the clumps (see 
Fig. |6j. In principle, rj determines to which extent the vorosity 
effect (i.e., the velocity gaps between the clumps, cf. lOwocki 
2008) is important for the line formation. As defined, rj does not 
contain any assumptions on the spatial structure of the wind. 
77 << 1 means that the velocity gaps between the clumps are 
much smaller than the thermal velocity, which in turn means that 
the probability for a photon to encounter a clump within its res- 
onance zone is high. If we assume each clump to be optically 
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thick, every encounter will lead to an absorption. Thus the prob- 
ability for photon absorption is high when the value of 77 is low. 
Vice versa, 77 >> 1 results in a high probability for the photon 
to escape its resonance zone without interacting with the wind 
material, i.e., a low absorption probability. If the entire velocity 
space were covered by clumps, 77 = 0. 

For the wind geometry used in our stochastic models, we 
may write (see AppendixlBlfor a derivation) 



vpSt(l - fy\Sv/Sv p \) St(l - f v \Sv/6v p \) v p 



v t 



,-2' 



(7) 



where L r is the radial Sobolev length of a smooth model, which 
for/3 = 1 is L r * v t r 2 (as usual, r and in /?* and St in fd yn )- Note 
that in Eq.Qalso the density-clumping parameters have entered 
the expression for 77, illustrating that there is an intimate coupling 
with the spatial clumping parameters, even though the vorosity 
effect initially depends on velocity parameters alone. For exam- 
ple, consider a wind with clumps that follow a smooth /3 velocity 
law. By bringing the clumps spatially closer together (for exam- 
ple by decreasing St), the velocity gaps between them decrease 
as well. Thus one may choose to describe the changed situa- 
tion either in terms of a less efficient porosity, because of fewer 
'density holes' in the reson ance zone through w hich the pho- 
tons can escape (as done by Osk inova etal]l2007h . or in terms 
of a less efficient vorosity, because of smaller velocity gaps be- 
tween the clumps. Of course, one may also obtain a lower ve- 
locity gap between the clumps by increasing the actual velocity 
spans inside the clumps, as simulated in our stochastic models 
when \dv/dvp\ > 1. This effect, leading to a rather low vorosity, 
has already been demonstrated to be at work in the RH models 
(Sect @3). 

Using the parameters of our default model, Fig. [6] (right 
panel) displays 77 as a function of velocity and shows that 77 in- 
creases rapidly in the inner wind, reaches a maximum at v « 
0.33, and then drops in the outer wind. To compare this behav- 
ior with that of the line profiles, we can associate absorption at 
some frequency jt b s with the corresponding value of the veloc- 
ity, because absorption occurs at x t, s ~ fiv ~ v (radial pho- 
tons dominate). In the default model's absorption-part line pro- 
file (see Fig. [4] the middle plot in the lower-left panel), a strong 
de-saturation occurs directly after the clumping is set to start (at 
r = 1.3, v « 0.23), followed by a maximum at x t, s ~ 0.35, and 
finally an absorption dip toward the blue edge. The behavior of 
the line profile is thus well mapped by 77, and we may explain 
the absorption dip as a consequence of the low value of 77 in the 
outer wind, which in turn stems from the slow variation of the 
velocity field (i.e., from radially extended resonance zones). 

5.2. Density parameters 

To isolate density-clumping effects, we use a smooth ft — 1 ve- 
locity law in this subsection. Despite the smooth velocity field, 
there are still holes in velocity space (because of the density 
clumping, at the locations where the ICM is present), and the 
expression for 77 (Eq. |7]i remains valid. Since a smooth velocity 
field corresponds to Sv - Svp, also the run of 77 is equal to the 
one displayed in Fig. [6] In this subsection we work only with in- 
tegrated profile strengths (characterized by the equivalent width 
Wa of the line's absorption part). The shapes of the line profiles 
are discussed in Sect. 16.11 

Fig. 13 shows Wa as a function of kq, for smooth models 
as well as for stochastic models with and without a contribut- 
ing ICM. The figure directly tells: i) The default model (x; c = 




0.01 0.10 1.00 10.00 100.00 

Fig. 7. Equivalent widths Wa of the absorption parts of line pro- 
files, normalized to the value of a saturated line, as a function 
of line strength parameter kq. The solid line is calculated from 
smooth models, and the dashed, dashed-dotted, and dotted lines 
from stochastic models with a smooth velocity field and St = 0.5, 
/ v = 0.25, and x K as indicated. 



0.0025) for the intermediate line {kq = 5.0) displays a Wa corre- 
sponding to a smooth model with a ko roughly ten times lower, 
ii) Lines never saturate if the ICM is (almost) void, iii) The run 
of Wa for the smooth and clumped models decouple well before 
kq reaches unity, iv) For intermediate lines, the response of Wa 
on variations of ko is weak for clumped models. Points one to 
three confirm our findings from Sect. 14.31 

A variation of 5t in the stochastic models affects primarily 
the high kq part (kq ^ 1.0) of the curves in Fig. [7] For example, 
lowering 5t in the model with a void ICM results in an upward 
shift of the dashed curve and vice versa. To obtain saturation 
with a void ICM, St * 0.05 is required, which may be under- 
stood in terms of Eq.|7] For St = 0.05, the 77-values correspond- 
ing to the default model are decreased by a factor of ten, and 77 
reaches a maximum of only about unity, with even lower values 
for the majority of the velocity space (cf. Fig. [6] right panel). 
The velocity gaps between the clumps then become closed, and 
the line saturates. In this situation, however, the intermediate line 
becomes saturated as well, again demonstrating the necessity of 
a non-void ICM to simultaneously saturate a strong line and not 
saturate an intermediate line. Only a properly chosen x ic param- 
eter ensures that the velocity gaps between the clumps become 
filled by low-density material able to absorb at strong line opac- 
ities, but not (or only marginally) at opacities corresponding to 
intermediate lines. 

When varying jq c , the primary change occurs at the high ko 
end of Fig. [7] For higher (lower) values of x; c , this part becomes 
shifted to the left (right), and the curve decouples earlier (later) 
from the corresponding curve for the void ICM. A higher ICM 
density obviously means that the ICM starts absorbing photons 
at lower line strengths and vice versa. Thus, observed saturated 
lines could potentially be used to derive the ICM density (or at 
least to infer a lower limit), if the mass-loss rate (and abundance) 
is known from other diagnostics. 

The behavior of the absorption with respect to the volume 
filling factor is as expected from the expression for 77; the higher 
/ v , the lower the value of 77, and the stronger the absorption. This 
is because a higher / v for a fixed St implies that the clumps be- 
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come more extended, whereas the distances between clump cen- 
ters remain unaffected. Consequently, a larger fraction of the to- 
tal wind velocity is covered by the clumps, leading to stronger 
absorption. For weak lines (ko as 0.05), the ratio W^/W^, sm devi- 
ates significantly from unity only when f v $0.1. Only for such 
low values can high enough clump densities be produced so that 
the clumps start to become optically thick. 

From Fig. [7J it is obvious that, generally, clumped models 
have a different (slower) response in W\ to an increase in Kq than 
do smooth models. This behavi or may be observa tionally tested 
using UV resonance doublets dMassa et al.l I2008T) . because the 
only parameter that differs between the two line components is 
the oscillator strength. Thus, if a smooth wind model is used 
and the fitted ratio of line strengths (i.e., /fo,t>iue/*'o,re<i) does not 
correspond to the expected ratio of oscillator strengths, one may 
interpret this as a signature of a clumped wind. Such behavior 
was found by Mas sa et al.l (|2008), where the observed ratios of 
the blue to red component of Si IV /i/il394,1403 inB supergiants 
showed a wide spread between unity and the expected factor of 
two. This result indicates precisely the slow response to an in- 
crease in kq that is consistent with inhomogeneous wind models 
such as those presented here, but not with smooth ones. In inho- 
mogeneous models, the expected profile strength (or Wa) ratio 
between two doublet components will depend on the adopted 
clumping parameters (as demonstrated by Fig.[7jand the discus- 
sion abov e) and ma y in principle take any value in the range 
found bv iMassa et all That is, while a profile-strength ratio de- 
viating from the value expected by smooth models might be a 
clear indication of a clumped wind, the opposite is not necessar- 
ily an indication of a smooth wind. Furthermore, the degeneracy 
between a variation of clumping parameters and kq suggests that 
un-saturated resonance lines should be used primarily as con- 
sistency tests for mass-loss rates derived from other diagnostics 
rather than as direct mass-loss estimators. We will return to this 
problem in Sect. 16.61 where a first comparison to observations is 
performed for the PV doublet. 

5.3. Velocity parameters 

The jump velocity parameter, Vj, affects only the strong lines (or, 
more specifically, the lines for which the ICM is significant), 
and determines the maximum velocity at which absorption can 
occur. For example, by setting vj = 0, no absorption at frequen- 
cies higher than x = 1 is possible (unless 5v is positive and 
very high). A higher Vj also implies more velocity overlaps, and 
thereby an increased amount of backscattering due to multiple 
resonance zones. Both effects are illustrated in Fig. [8] Judging 
from the line profiles of the lower panel, the blue edge and the 
reduction of the re-emitted flux blueward of the line center may 
both be used to constrain vj . The upper panel shows one slice of 
the corresponding velocity fields, illustrating that the underlying 
/3 law is recovered almost perfectly when using Vj = O.Olv^g and 
5v - dvp. With this velocity law and a non-void ICM, the cor- 
responding strong line profile is equivalent to a profile from a 
smooth model. 

In Sects. 1431 and 14.41 we showed that a higher value of the 
clumps' velocity spans led to stronger absorption for interme- 
diate lines. In principle this is as expected from Eq. [7] where 
77 always decreases with increasing \8v/Svp\. However, with the 
very high value of \6v/6vg\ used in, e.g., the RHcopy model, 
one realizes that 77 in Eq. |7j becomes identically zero, because 
fv\6v/5vp\ - 1. An 77 — corresponds to the whole velocity 
space being covered by clumps, and the saturation limit should 
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Fig. 8. Upper: Velocity structures (one slice) in two stochas- 
tic models with density-clumping parameters as for the default 
model, and different velocity parameters. Dashed: dv/dvp = 1 
and Vj/Vjg = 0.01. Solid: Sv/dvp = — 1 and vjvp = 0.5 below 
vg = 0.6 and Vj/v/j = 0.15 above. Lower: Corresponding line 
profiles for a strong line. 



be reached. As is clear from Fig. |4] however, this is not the case. 
This points out two important details not included when deriv- 
ing the expression for 77 and interpreting the absorption in terms 
of this quantity, namely that clumps are distributed randomly 
(with 6t determining only the average distances between them) 
and that the parameter Vj allows for an asymmetry in the ve- 
locities of the clumps' starting points (see Sect. I2.21 >. These two 
issues lead to overlapping velocity spans for some of the clumps, 
whereas for others there is still a velocity gap left between them, 
through which the radiation can escape. Therefore the profiles 
do not reach complete saturation, despite that on average 77 = 0. 
This illustrates some inherent limitations when trying to inter- 
pret line formation in terms of a simplified quantity such as 77. 

The impact from the velocity spans of the clumps on the 
line profiles also depends on the density-clumping parameters. 
To achieve approximately the same level of absorption, a higher 
value of 5v/5vp was required in the RHcopy model (/ v = 0.1) 
than in the default model (/ v = 0.25), see Fig.|4] Since Svp °c f v 6t 
(see AppendixlBl. the actual velocity spans of the clumps are dif- 
ferent for different density-clumping parameters, even if 6v/Svp 
remains unchanged. 

By changing the sign of 6v in the default model (that is, as- 
suming a positive velocity gradient inside the clumps), we have 
found that our results qualitatively depend only on |<5v|. Some 
details differ though. For example, a Sv > in our stochastic 
models permits absorption at velocities higher than the terminal 
one also within the clumps, whereas Sv < restricts the clump 
velocities to below the local vp (see Fig . [2J . In this matter vj plays 
a role as well, since Vj controls where, with respect to the local 
Vfi, the clumps begin. For reasonable values of Vj, however, its in- 
fluence is minor on lines where the ICM is insignificant. Finally, 
tests have confirmed that optically thin lines are only marginally 
affected when varying Sv/Svp. 
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Fig. 9. Total, absorption part, and re-emission part line profiles 
for ID, smooth models with = 5.0 (dashed-dotted lines) and 
kq = 5.0/(277) (solid lines, see Sect. 16. U . and for a 2D, stochastic 
model with density parameters as the default model and a /3 = 1 
velocity law (dashed lines). 
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Fig. 10. Upper panel: Density structures of one slice in the de- 
fault stochastic model (upper), in the default stochastic model 
with a modified St (middle, see Sect. 16.21 1. and in FPP (lower). 
Lower panel: Line profiles for the absorption part of an inter- 
mediate line for the default model (solid line), for the default 
model with a modified 5t (dashed line), and for the default model 
with an ionization structure decreasing with increasing velocity 
(dashed-dotted line, see text). 



6. Discussion 

6.1. The shapes of the intermediate lines 

For intermediate lines, the shape of the absorption part of the 
default model differs significantly from the shape of a smooth 
model (see Fig. |4] the middle plot in the lower-left panel). We 
showed in Sect. I5.ll that the shapes could be qualitatively un- 
derstood by the behavior of 77. This is further demonstrated here 
by scaling the line strength parameter of a ID, smooth model, 
using a parameterization /cq oc rf x outside the radius r — 1.3 
where clumping is assumed to start. Fig.|9]displays the line pro- 
files of ID, smooth models with Kq = 5.0 and Kq = 5.0/(2rf). 
These profiles are compared to those calculated from a 'real' 2D 
stochastic model with density-clumping parameters as the de- 
fault model, but with a /3 = 1 velocity field. 77 was calculated 
from Eq. [7] using the parameters of the default model and a 
/3 = 1 velocity law, and the factor of 2 in the denominator of 
the scaled kq was chosen so that the integrated profile strength 
of the 2D model was roughly reproduced. From Fig.|9]it is clear 
that the ID model with scaled kq well reproduces the 2D results, 
indicating that indeed 77 governs the shape of the line profile. We 
notice also that these profiles display a completely black absorp- 
tion dip in the outermost wind, as opposed to the default model 
with a non-monotonic velocity field (see Fig. |4] the middle plot 
in the lower-left panel). This is because the /3 velocity field does 
not allow for any clumps to overlap in velocity space (see the 
discussion in Sect. 15.3b . making the mapping of 77 almost per- 
fect. 

Let us also point out that the line shapes can be somewhat 
altered by using a different velocity law, e.g., /3 + 1. Such a 
change would affect the distances between clumps as well as the 
Sobolev length, and thereby the line shapes of both absorption 
and re-emission profiles. However, in all cases is the shape of 
the re-emission part similar in the clumped and smooth models. 



6.2. The onset of clumping and the blue edge absorption dip 

We have used r = 1.3 as the onset of wind clumping in 
our stochastic models, which roughly corresponds to the radius 
where significant structure has developed from the line - driven 
instab ility in our RH models. However, iBouret et al ] (I2003L 
2005) analyzed O-stars in the Galaxy and the SMC, assum- 
ing optically thin clumps, and found that clum ping star t s deep 
in the wind, just above the sonic point. Also iPuls et all {2006) 
used the optically thin clumping approach, on p 2 -diagnostics, 
and found similar results, at least for O-stars with dense winds. 
With respect to our stochastic models, the qualitative results 
from Sects. [4]and[5]remain valid when choosing an earlier onset 
of clumping. Quantitatively, the integrated absorption in inter- 
mediate lines becomes somewhat weaker, because the clumping 
now starts at lower velocities, and of course the line shapes in 
this region are affected as well. The onset of wind clumping will 
be important when comparing to observations, as discussed in 
Sect. EH 

The stochastic models that de-saturate an intermediate line 
generally display an absorption dip toward the blue edge (see 
Figs.[4]and[9]i, which has been interpreted in terms of low values 
of 77 in the outer wind (see Sect. 15. 11 1. However, this characteris- 
tic feature (not to be confused with the so-called DACs, discrete 
absorption components) is generally not observed, and one may 
ask whether it might be an artifact of our modeling technique. 
In the following we discuss two possibilities that may cause our 
models to overestimate the absorption in the outer wind; the ion- 
ization fraction and too low clump separations. 

Starting with the former, we have so far assumed a constant 
ionization factor, q — 1 (cf. Eq. IA.9b . This is obviously an over- 
simplification. For example, an outwards decreasing q would re- 
sult in less absorption toward the blue edge. Here we merely 
demonstrate this general effect, parameterizing q = vo/vp in the 
stochastic default model (see Table with vo = 0.1 the starting 
point below which q = 1. Fig. [TUl (lower panel, dashed-dotted 



J.O. Sundqvist et al.: Mass loss from inhomogeneous hot star winds 



11 



lines) shows how the absorption in the outer wind becomes sig- 
nificantly reduced. 

The temperature structure of the wind is obviously important 
for the ionization balance. Whereas an isothermal wind is as- 
sumed in POF (see Sect. 12. Il l, the FPP model has shocked wind 
regions with temperatures of several million Kelvin. To roughly 
map corresponding effects on the line profiles, we re-calculated 
profiles based on FPP models assuming q — in all regions with 
temperatures higher than T — 10 5 K, and q — 1 elsewhere. Since 
the hot gas resides primarily in the low-density regions, how- 
ever, the emergent profiles were barely affected, and particularly 
intermediate lines remained unchanged. 

On the other hand, the X-ray emission from hot stars (be- 
lieved to originate in clump-clump collisions, see FPP) is known 
to be crucial for the ionization balance of highly ioniz ed species 
such as C IV, N V, and O VI (see, e.g., the discussion in Puis et al] 
2008b). X-rays have not been included here, but could in prin- 
ciple have an impact on our line profiles, by illuminating the 
over-dense regions and thereby changing the ionization balance. 
Krticka&Kubaddl009), however, find that incorporating X-rays 
does not influence the PV ionization significantly. Finally, non- 
LTE analyses including feedback from optically thin clumping 
have shown that this as well has a s ignificant effect on t he derived 
ioniza tion fractions of, e.g., PV (Bou ret et alJl2005t iPuls et ail 
2008a). To summarize, it is clear that a full analysis of ionization 
fractions must await a future non-LTE application that includes 
relevant feedback effects from an inhomogeneous wind on the 
occupation numbers. 

In RH models, the average distance between clumps in- 
creases in the oute r wind, due to clump-clump collisions and 
veloc ity stretching (Feldmeie r et al.l[l997t [kunacres &~Owockil 
2002). Neglecting the former effect, our stochastic models have 
clumps much more closely spaced in the outer winc0. We have 
therefore modified the default model by setting 5t — 3 outside a 
radius corresponding to vp = 0.7. This is illustrated in the upper 
panel of Fig. [10] The mass loss in the new stochastic model is 
preserved (because the clumps are more extended, see the fig- 
ure), and this model now better resembles FPP. Recall that dif- 
ferences in the widths of the clumps are expected, since in the 
default model / c i w f' 1 = 4, whereas in FPP f c \ « 10. The corre- 
sponding line profile shows how the absorption outside x ~ 0.7 
has been reduced, as expected from the higher St. 

6.3. The velocity spans of the clumps 

In Sect. 14.41 it was found that \Sv\ > Svp in the RH models. 
Fig. [TTJ upper panel, shows the velocity spans of density en- 
hancements (identified as having a density higher than the cor- 
responding smooth value) in the FPP model, and demonstrates 
that, after structure has developed, \6v\ is much higher than Svp 
throughout the whole wind. These high values essentially stem 
from the location of the starting points of the density enhance- 
ments, which generally lie before the velocities have reached 
their post shock values (see Fig. QT| middle and lower panels). 
By using a /3 velocity law (which in principle corresponds to a 
stochastic velocity law with Vj = and 5v = Svp, see Fig. [8} 
together with the density structure from FPP, we simulated a RH 
wind with low velocity spans. Indeed, for the corresponding in- 
termediate line the equivalent width of the absorption part was 
~ 35% lower than that of the original FPP model. The strong 
line, on the other hand, remained saturated, because the ICM in 



5 The effect is minor in POF, since these RH models only extend to 
r~ 5 (see Sect. EH. 
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Fig. 11. Upper: Velocity spans of density enhancements in the 
FPP model (squares) and corresponding f3 intervals (diamonds). 
Lower: Three density enhancements and corresponding velocity 
spans in the FPP model, highlighted as in Fig.Q] 

FPP is not void. So, again, the RH models would in parallel dis- 
play de-saturated intermediate lines and saturated strong lines, 
were it not for the large velocity spans inside the clumps. 

We suggest that the large velocity span inside a shell (clump) 
is primarily of kinematic origin, and reflects the formation his- 
tory of the shell. The shell propagates o utwards through th e 
wind, essentially with a/3 — 1 velocity law (Owo cki et alj! 988). 
Fast gas is decelerated in a strong reverse shock at the inner rim 
of the shell. The shell collects ever faster material on its way out 
through the wind. This new material collected at higher speeds 
resides on the star-facing side, i.e. at smaller radii, of the slower 
material collected before. Thus, a negative velocity gradient de- 
velops inside the shell. The fact that |<5v| » Svp in FPP seems 
to reflect that the shell is formed at small radii, and then advects 
outwards maintaining its steep interior velocity gradient- From 
this formation in the inner, steeply accelerating wind, velocity 
spans within the shells up to (a few) hundred kms , as seen in 
Fig-Qj] appear reasonable. 

However, the dynamics of shell formation in hot star winds 
is very complex due to the creation and subsequent merging of 
subshells, as caused by nonlinear perturbation growth and the re- 
lated excitation of ha rmonic overtones of the perturbation period 
at the wind base (see iFeldmeierl [19951) . Future work is certainly 
needed to clarify to which extent the large velocity spans inside 
the shells in RH models are a stable feature (see also Sect. 17. 2\ . 

6.4. 3D effects 

A shortcoming of our analysis is the assumed symmetry in <t>. 
The 2D rather than 3D treatment has in part been motivated 
by computational reasons (see Appendix [All. More importantly 
though, we do not expect our qualitative results to be strongly 
affected by an extension to 3D. Within the broken-shell wind 
model, all wind slices are treated independently, and distances 
between clumps increase only in the radial direction. Therefore 



6 Actually, the velocity gradient may further steepen during advec- 
tion, due to faster gas trying to overtake slower gas ahead of it; however, 
this effect is balanced by pressure forces in the subsonic postshock do- 
main. 
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the expected outcome from extending to 3D is a smoothing effect 
rather than a reduction or increase in integrated profile strength 
(similar to the smoothing introduced by N@, see Sect. 14. U . Also, 
we have shown that the main effect from the inhomogeneous 
winds is on the absorption part of the line profiles (see, e.g., 
Sect. 16. U . The formation of this part is dominated by radial pho- 
tons, especially in the outer wind, because of the dependence 
only on photons released directly from the photosphere. This 
implies that most photons stay within their wind slice, restricting 
the influence from any additional 'holes' introduced by a broken 
symmetry in <1> to the inner wind. Of course, these expectations 
hold only within the broken shell model, because in a real 3D 
wind the clumps will, for example, have velocity components 
also in the tangential directions. 

6.5. Comparison to other studies 

To scale the smooth opaci ty in the formal integra l of the non-LTE 
atmospheric code Po WR, lOskinova et al.l d2007) used a porosity 
formalism in which both / v and the average distance between 
clumps enter. Other assumptions were a void ICM, a smooth 
ft velocity field, and a microturbulent velocity v t ss 50kms _1 , 
the last identified as the velocity dispersion within a clump. 
However, a direct comparison between their study and ours is 
hampered by the different formalisms used for the spacing of the 
clumps. Here we have used the 'broken-shell' wind model as a 
base (see Sect. 12.21 1. in which each wind slice is treated indepen- 
dently and the distance between clumps increases only in the ra- 
dial direction (clumps preserve their lateral angles). This gives a 
radial number density of cl umps, n z \ oc v , the same as used by, 



e.g.,[(Dskinova et al. (20 0^), when synth esizing X-ray emission 
from hot stars. In Oskinova et al. (2007), on the other hand, the 
distance between clumps increases in all spatial directions. In a 
spherical expansion, this gives a radial number density of clumps 
n d oc v~'r~ 2 , i.e., clumps are distributed much more sparsely 
within this model, especially in the outer wind. Therefore their 
choice of Lq = 0.2 is not directly comparable with 6t = 0.2 in our 
models. The shapes of the clu mps differ between the two mod- 
els as well; in lOskinova et al.l clumps are assumed to be 'cubes', 
whereas here the exact shapes of the clumps are determined by 
the values of the clumping parameters. Despi te these di fferences, 
our findings confirm the qualitative results of lOskinova et alJ that 
the line profiles become weaker with an increasing distance be- 
tween clumps as well as with a decreasing v t . These results may 
be interpreted on the basis of the effective escape ratio, rj (see 
Eq. [7j. Both a decrease in \\ and an increase in the distance be- 
tween clumps mean that the velocity span covered by a reso- 
nance zone becomes smaller when compared to the velocity gap 
between two clumps (see Fig. [6] left panel), leading to higher 
probabilities for line photons to escape their resonance zones 
without interacting with the wind material. 

An important result of this paper is that models that de- 
saturate intermediate lines require a non-void ICM t o saturate 
strong lines. This is confirmed by the Oskinova et al.l model, in 
which the ICM is v oid and strong lines indeed do not saturate 
dHamann et al.ll2009t) . 

lOwockil (120081) proposed a simplified description of the non- 
monotonic velocity field to account for vorosity, i.e., the veloc- 
ity gaps between the clumps. Here, the vorosity effect has been 
discussed using the quantity rj (see Sect. I5.lt . and we have in- 
troduced two new parameters to characterize a non-monotonic 
velocity field, 6v and Vj. The reason for introducing a new pa- 
rameterization is that when using a single velocity parameter, we 
have not been able to simultaneously meet the constraints from 



strong, intermediate, and weak lines as listed in Sect. |4] Tests 
using a ' veloci ty clumping factor' / ve i = Sv/Av as proposed by 
lOwockil (|2008), together with a smooth density structure, have 
shown that this treatment indeed can reduce the line strengths 
of intermediate lines, but that the observational constraints from 
strong lines may not be met. Still, the basic concept of vorosity 
holds within our analysis. For example, one may phrase the high 
values of 6v in the RH models in terms of insufficient vorosity. 

6.6. Comparison to observations 

We finalize our discussion by performing a first comparison to 
observations. The two components of the Phosphorus V AAl 118- 
1 128 doublet are rather well separated, and the singlet treatment 
used here suffices to model the major part of the line complex. 
Nevertheless, the two components overlap within a certain re- 
gion (indicated in Fig. IT2b . so when interpreting the results of 
this subsection, one should bear in mind that the overlap is not 
properly accounted for, but treated as a simple multiplication of 
the two profiles. 

We used observed FUSE spectra (kindly provided by A. 
Fullerton) from HD 210839 (A Cep), a supergiant of spectral 
type 06 1(n)fp. When computing synthetic spectra, we first as- 
sumed optically thin clumping with a constant clumping fac- 
tor / c i = 9 and a smooth f3 = 1 velocity field, fg = 9 agrees 
fairly well with the analysis of iPuls et al.l d2006l) . who derived 
clumping factors f c \ - 6.5 for r ~ 1.2... 4.0 and f c \ = 10 
for r ~ 4.0 ... 15, assuming an un-clumped outerm ost windf! 
We too k the ionization fraction q = q(r) of PV from Pui s et al.l 
(l2008al) . calculated with the unified non-LTE atmosphere code 
FASTWIN D for an Q6 supergian t, using the Phosphorus model 
atom from Pauldra ch et al.l ([2001). The feedback from optically 
thin clumping was accounted for and X-rays were neglected. 
This ionization fraction was then used as input in our MC- 
1D code when computing the synthetic spectra. We assigned 
a thermal plus a highly supersonic 'microturbulent' velocity 
Vt = 0.05 (corresponding to 110 kms -1 ), as is conventional in 
this approach. The mass-loss rate was derived using the well 
known relation between /<o and M (e.g.. |Puls et al.l l2008b). For 
atomic and stellar pa rameters, we adopted the same values as in 
iFullerton et all(l2006l) . 

The dashed line in Fig. [T2l represents our fit to the observed 
spectrum, assuming optically thin clumping, resu lting in a mass- 
loss ra te M = 0.24, in units of 10~ 6 M^ vr~' . [Fullerton et al.l 
(2006) derived {q)M — 0.23 for this star. Because our clumped 
FASTWIND model predicts an averaged ioni zation fraction 
(q) ~ 0.9 in the velocity regions utilized by Full erton et all 
the two rates a re in excellent agreement. On the other hand, 
iRepolust et all (120041) for HD 210839 derived M = 6.9 from 
H a assuming an unclumped wind, yielding Mh„ = 2.3 when 
accounting for the reduction implied by our assumed f c \ - 9 

— 1/2 

(Muc, - Mu a ,smf cl )• This rate is almost ten times higher than 
that inferred from PV, and thus results in PV line profiles that 
are much too strong (see Fig.[T2] dashed-dotted line). That is, to 
reconcile the H a and PV rates for HD 210839 with models that 
assume optically thin clumps also in PV, we would have to raise 
the clumping factor to f c \ > 100. In addition to this very high 
clumping factor, the low rate inferred from the PV lines con- 
flicts with the theoretical value M = 3.2 provided by the mass- 



7 This stratification has been found to be prototypical for O- 
supergiants and was, together with its well developed PV P Cygni pro- 
files, the major reason for choosing A Cep as comparison object instead 
of, e.g., f Pup, which displays a somewhat unusual run of f c \. 
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loss recipe in Vink et afl (120001) (using the stellar parameters of 
iRepolust et all 20041) . and is also strongly disfavored by current 
massive star evolutionary models (lHirschill2008l) . 

Next we modeled the PV lines using our MC-2D code to- 
gether with a stochastic 2D wind model. The same clump- 
ing factor (/ c i = 9) and ionization fraction (calculated from 
FASTWIND, see above) were used. This time, we assigned 
v t = 0.005, i.e., applied no microturbulence. In previous sec- 
tions, e.g. 14.31 and 16.11 we showed that stochastic models gen- 
erally display a line shape different from smooth models, with 
a characteristic absorption dip at the blue edge as well as a dip 
close to the line center. Such shapes are not seen in the PV lines 
in A Cep. Thus, to better resemble the observed line shapes, we 
used different values for St and x; c in the inner and outer wind 
(the former modification already discussed in Sect. 16.21 ) and let 
clumping start close to the wind base. Clumping parameters are 
given in Table|2] model Obsl. 

As illustrated in Fig. Q~2] the synthetic line profiles using 
M = 2.3, as inferred from H a , are now at the observed levels. 
Because of our insufficient treatment of line overlap, we gave 
higher weight to the Al 1 18 component when performing the fit- 
ting, but the profile-strength ratio between the blue and red com- 
ponent was nevertheless reasonably well reproduced (see also 
discussion in Sect. 15.21 ). However, though the fit appears quite 
good, we did not aim for a perfect one, and must remember the 
deficits of our modeling technique. For example, while the early 
onset of clumping definitely improved the fit (using our default 
value, there was a dip close to line center) and might be consid- 
ered as additional evidence that clumping starts close to the wind 
base, the same effect could in principle be produced by non-LTE 
effects close to the photosphere or by varying the underlying /3 
velocity law. Such effects will be thoroughly investigated in a 
follow-up paper, which will also include a comparison to obser- 
vations from many more objects. 

Clearly, a consistent modeling of resonance lines (at least 
of intermediate strengths) requires the consideration of a much 
larger parameter set than if modeling via the standard diagnos- 
tics assuming optically thin clumping, and a reasonable fit to a 
single observed line complex can be obtained using a variety 
of different parameter combinations. The analysis of PV lines 
as done here can therefore, at present, only be considered as 
a consistency check for mass-loss rates derived from other, in- 
dependent diagnostics, and not as a tool for directly estimating 
mass-loss rates. Additional insight might be gained by exploiting 
more resonance doublets, due to the different reactions of profile 
strengths and shapes on ko. The different slopes of the equiva- 
lent width as a function of /cq in smooth and clumped models, 
especially at intermediate line strengths (Sect. [52]), may turn out 
to be decisive. However, because of, e.g., the additional impact 
from the ICM density, also this diagnostics requires additional 
information from saturated lines. Taken together, only a consis- 
tent analysis using different diagnostics and wavelength bands, 
and embedded in a suitable non-LTE environment, will (hope- 
fully) provide a unique view. 

7. Summary and future work 

7.1. Summary 

Below we summarize our most important findings: 

- When synthesizing resonance lines in inhomogeneous hot 
star winds, the detailed density structure, the non-monotonic 
velocity field, and the inter-clump medium are all important 
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Fig. 12. Observed FUSE spec tra of the PV double t AAl 1 18-1 128 
for the 06 supergiant A Cep (iFullerton et al.l2 006). The synthetic 
spectra are calculated for two ID models assuming optically thin 
clumping (see Sect. 16.61 ) and for one 2D stochastic model with 
parameters as in Tabled model Obsl. The models have mass- 
loss rates M [M yr -1 ] as given in the figure. The zero point fre- 
quency is shifted to the line center of the Al 1 18 component, and 
the two arrows at the bottom of the figure indicate in which re- 
gion the two components overlap. 



for the line formation. Adequate models must be able to si- 
multaneously meet observational and theoretical constraints 
from strong, intermediate, and weak lines. 
Resonance lines are basically unaffected by the inhomoge- 
neous wind structure in the limit of optically thin clumps, but 
the clumps remain optically thin only for very weak lines. 
We confirm the basic effects of porosity (stemming from op- 
tically thick clumps) and vorosity (stemming from velocity 
gaps between the clumps) in the formation of primarily lines 
of intermediate strengths. 

We point out the importance of a non-void ICM for the si- 
multaneous formation of strong and intermediate lines that 
meet observational constraints. 

Porosity and vorosity are found to be intrinsically coupled 
and of similar importance. To characterize their mutual effect 
on intermediate lines, we have identified a crucial parameter, 
the 'effective escape ratio', that describes to which extent 
photons may escape their resonance zones without ever in- 
teracting with the wind material. 

We confirm previous results that time-dependent, radiation- 
hydrodynamic wind models reproduce observed characteris- 
tics for strong lines, without applying the highly supersonic 
microturbulence needed in smooth models. 
A significant profile strength reduction of intermediate 
lines (as compared to smooth models) is for the radiation- 
hydrodynamic models prevented by the large velocity spans 
of the density enhancements, implying that the wind struc- 
tures predicted by present day RH models are not able to re- 
produce the observed strengths of intermediate lines unless 
invoking a very low mass-loss rate. 

Provided a non-void ICM and not too large velocity spans in- 
side the clumps, 2D stochastic wind models saturate strong 
lines, while simultaneously not saturating intermediate lines 
(that are saturated in smooth models). Using typical vol- 
ume filling factors, / v * 0.25, the resulting integrated profile 
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strength reductions imply that these inhomogeneous models 
would be compatible with mass loss rates roughly a factor 
of ten higher than those derived from resonance lines using 
smooth models. 

- A first comparison to observations was made for the 06 su- 
pergiant A Cep. It was found that, indeed, the line profiles 
of PV based on a 2D stochastic wind model, accounting for 
a detailed density structure and a non-monotonic velocity 
field, reproduced the observations with a mass-loss rate al- 
most ten times higher than the rate derived from the same 
lines, but with a model that used the optically thin clumping 
approach. This alleviated the discrepancies between theoreti- 
cal predictions, evolutionary constraints, and previous mass- 
loss rates based on winds assumed either to be smooth or to 
have optically thin clumps. 

7.2. Future work 

We have investigated general properties of resonance line for- 
mation in inhomogeneous 2D wind models with non-monotonic 
velocity fields. To perform a detailed and quantitative compar- 
ison to observations, and derive mass-loss rates, simplified ap- 
proaches need to be developed and incorporated into non-LTE 
models to obtain reliable occupation numbers. Extending our 
Monte-Carlo radiative transfer code to include line overlap ef- 
fects in doublets is critical for more quantitative applications, 
and an extension to 3D is also necessary. Further applications 
involve synthesizing emission lines, for example to test the opti- 
cally thin clumping limit both in the parameter range where this 
is thought to be appropriate (e.g., for 0-/early B-stars), and in 
other more complicated situations. Indeed, the present genera- 
tion of line-blanketed model atmospheres does not seem to be 
able to reproduce H a line profiles from A-supergiants, which 
are observed as P-Cygni profiles with non-saturated troughs, 
whereas the simulations (assuming optically thin clumping) re- 
sult in saturated troughs (R.-P. Kudritzki, private communica- 
tion). Since H,,, is a quasi-resonance line and not a recombina- 

i ii "~i 

tion line in these cooler winds (e.g., Kudritzki & Puis 2000), this 

behavior might be explained by the presence of optically thick 
clumps. 

Finally, it needs to be clarified if the large velocity span in- 
side clumps generated in RH models is independent of addi- 
tional physics that is not, or only approximately, accounted for 
in present simulations (such as more-D effects and/or various 
exciting mechanisms). If the large velocity span is a stable fea- 
ture, one might come to the (rather unfortunate) conclusion that 
either the observed clumping features are not, or only weakly, re- 
lated to the line-driven instability, or the discrepancies between 
observed and synthetic flux distribution (from the X-ray to the 
radio regime) might involve processes different from the present 
paradigm of wind clumping. 
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Appendix A: The Monte-Carlo transfer code 

A. 1 . The code 

Here we describe our Monte-Carlo radiative transfer code (MC- 
2D) in some detail. For an overview of basic assumptions, see 
Sect.|3]in the main paper. For testing purposes, versions to treat 
spherically symmetric winds, either in the Sobolev approxima- 
tion (MCS-1D) or exactly (MC-1D), have been developed as 
well. 

Geometry. For wind models in which the spherical symmetry is 
broken, we can no longer restrict photon trajectories to rays with 
constant impact parameters (see below). Moreover, the observed 
spectrum will depend on the observer's placement relative to the 
star. Fig. IA. 1 1 illustrates the geometry in use, a standard right- 
handed spherical system (r, 0, <E>) defined relative to a Cartesian 
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set (X, Y, Z) (transformations between the two may be found in 
any standard mathematical handbook). At each coordinate point 
we also construct a local coordinate system using the local unit 
vectors (r u , @ u , <L> U ), which for a photon propagating in direction 
n a is related to the radiation coordinates (9, <p) (see Fig. IA.lb via 



cos 6 = p. = r u 



sin d> sin 6 — O u • n u = 
cos sin 6 — U ■ n u = [O u X r u ] ■ n u . 



(A.l) 
(A.2) 
(A.3) 



The radiation coordinates are defined on the intervals 8 — . . . n 
and <p = . . . 2jt, but due to the symm etry in O, only the range 
(/> — . ..7i needs to be considered (see iBusche & Hillierl 2000). 
Also, for this symmetry, the direction cosines of n u simplify to 



fx sin + yjl - p 2 cos <p cos 0, 



"y = V^"~ 



/j. 2 sin ( 



n z - fi cos © - ^1 - 



p 2 COS ( 



sin 0. 



(A.4) 



(A.5) 



(A.6) 



Eqs. IA.1HA.6I are used to update the physical position (r, 0) 
of the photon and the local values of the radiation coordinates 
(9, (p). By tracking the photon on a radial mesh, both the physical 
and radiation coordinates can be updated exactly. Interpolations 
are necessary only when a photon is scattered or when it crosses 
a 0-boundary to another win d slice. Essentially the s ame coor- 
dinate system is used by, e.g.. lBusche & Hillierl d2000l) . We col- 
lect escaped photons according to their ©-angles at 'infinity^ 
and bin them using the same N@ bins as in the underlying wind 
model (see Sect. |2j. 

For spherically symmetric wind models, we adhere to the 
customary (p,z) spatial coordinate system with p being the im- 
pact parameter and z the direction toward the observer. Each time 
a photon is scattered and its direction determined, a new impact 
parameter is computed from the relation p — r V(l - p 2 ), ap- 
preciating that all points on a surface of constant radius can be 
treated equally in this geometry. 

Releasing photons. We release photons from the lower bound- 
ary uni formly in <j> and with a distribution function oc pdp in 
p. (e.g.. lLucvm 983). The angular coordinate is selected so 
that photons are uniformly distributed over the surface area 
dA = sin ®d®d<&. 



Absorption. The probability of photon absorption is oc e~ T dr, 
hence the optical depth r the photon travels before absorption 
can be selected according lo t - - h\R\, where R\ is a random 
number between and 1 . The position for absorption in the wind 
may then be determined by inverting the line optical depth inte- 
gral along the photon path 



t v = J Xvds, 

with the frequency-dependent opacity 



Xv = KLP<Pv 



(A.l) 



(A.8) 



8 The full 3D problem would require binning in O as well, which in 
turn would require a large increase in the number of simulated photons. 




Fig. A.l. Illustration of the coordinate system, see text. A color 
version of this figure is available in the web version. 



with <p v the absorption profile, kl the frequency integrated mass 
absorption coefficient, and p the mass density. All dependen- 
cies on spatial location are for simplicity suppressed here and 
in th e following. Fo r the opacity we use the parameterization 
from lHamannl (1 19811) and POF, 



K h Ap = -r Kopq, 

M 



(A.9) 



where A is the wavelength of the considered transition, kq is a 
'line-strength' parameter taken to be constant, M the radially and 
laterally averaged mass-loss rate, and q - q(r, 0) the fraction 
of the considered element that resides in the investigated ionic 
stage. Default here is q = 1, but effects from other ionization 
structures are discussed in Sect. 16.21 kq is proportional to the 
product of mass-loss rate and abundance of the considered ion, 
and, for a smooth wind, kq — 1 and kq = 100 give a typical 
medium and strong line, respectively. The parameterization as 
defined in Eq. |A.9| has the advantage that for smooth winds the 
radial optical depth in the Sobolev approximation collapses to 



T"Sob = 



r 2 vdv/dr ' 



(A. 10) 



when v and r are expressed in normalized units. The correspond- 
ing expression for clumpy winds is provided in Eq. The ab- 
sorption profile is assumed to be a Gaussian with a Doppler 
width Vt that contains the contributions from thermal and (if 
present) 'microturbulent' velocities. To solve Eq. IA.7l we adopt 
the dimensionless frequency x with the terminal velocity of a 
smooth outflow as the reference speed, 



v - Vp c 



(A. 11) 



and transform to the co-moving frame (hereafter CMF). vq is the 
rest-frame frequency of the line center and c the speed of light. 
We now assume that between two grid points the variation of the 
factor klpI\Q\ (see below) is small and may be replaced by an 
average value. The optical depth At v between two subsequent 
spatial points (r, 0) then becomes 



/LR* K L p -Aerf[x cmf /v t ] 
At v = | — x 1 

Voa Q 2 



(A. 12) 
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where Aerf is the difference of the error-function between the 
points, x cm f the dimensionless CMF frequency, and v, is calcu- 
lated in units of v x . Q = « u • V (n u ■ v) is the local directional 
derivative of the velocity in direction n u , with velocities mea- 
sured in units of and radii in units of By interpolating to 
the border whenever a photon crosses a © boundary, we locally 
recover the spherically symmetric expression 

<9v -> v i 
Q = -fi 2 + -(1 - /j 2 ). (A.13) 
or r 

For spherically symmetric winds, we have written a second im- 
plementation that allows for line transfer using the Sobolev ap- 
proximation. With this method each resonance zone is approx- 
imated by a point and the line only collects optical depth at at- 
mospheric locations where the observer's frame frequency x b s 
has been Doppler shifted to coincide with the CMF frequency 
for the line center. The condition for interaction thus is x b s = fiv 
and the last factor in Eq. lA.12l collapses to unity when calculat- 
ing the Sobolev optical depth. The Sobolev approach can be ex- 
pected a reasonable approximation when the variation of the fac- 
tor klp/\Q\ is small within the whole resonance zone contribut- 
ing to the optical depth in Eq. lA.12| i.e., small on length scales at 
least a few times the Sobolev length L = v t /\Q\. However, also in 
the Sobolev approximation more than one resonance point may 
be identified in a wind with a non-monotonic velocity field. 
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Fig. A.2. Synthetic line profiles for spherically symmetric mod- 
els, calculated with the labeled methods. Profiles are shown for 
a smooth model with acq = 1.0 and v t = 0.2 (upper) and for two 
POF snapshots with kq = 100 (middle) and kq =5.0 (lower) and 
v t = 0.005. The 2D profile is for an observer at the equator, x is 
the normalized observer's frame frequency (see Eq. IA.1U , and 
the ordinate displays the emergent flux normalized to the contin- 
uum flux. 



Re-emission. We assume complete redistribution and isotropic 
re-emission in the CMF, allowing for a multitude of scatter- 
ing events within one resonance zone. When the Sobolev ap- 
proximation is applied, re-emission is assumed to be coherent 
in the CMF and for the angular re- distribution w e then use the 
corresponding escape probabilities (ICastorll 19701). corrected for 
a treatment of negative velocity gradients (Rvb icki & Hummer! 
1978; POF). In this case, there is only one effective scattering 
event inside the localized resonance zone. 

After the photon has been re-emitted at some atmospheric 
location, the procedure runs again and searches for another ab- 
sorption. 

A.2. Radiative transfer code tests 

In this subsection we describe some of the verification tests of 
our MC radiative transfer code that we have made. The MC- 
ID version was first applied on spherically symmetric winds, 
comparing profiles from smooth, stationary win ds to profiles 
calcul ated usin g the well-established CMF (cf. iMihalas et al.l 
1975: lHamannl 1981) and SEI methods, and profiles from time- 
dependent RH winds to profiles calculated using the Sobolev 
method developed in POF. Thereafter we applied the MC-2D 
version on models in which all lateral slices had the same radial 
structure, comparing the results to the MC-ID version. 

First we calculated line profiles for smooth, ID winds. We 
have verified that for low0 values of v t , profiles from all the meth- 
ods described above agree perfectly, whereas for higher values 
the MC-ID and CMF give identical results but the SEI deviates 
significantly, especially for a medium-strong line (see Fig. IA.2l 
upper panel). This is due to the hybrid nature of the SEI tech- 
nique, which approximates the source function with its local 
Sobolev value but carries out the exact formal integral. Because 
of this, the method does not account for the increasing amount of 

9 For a typical terminal velocity value v„ = 2000 km s _1 , v, = 0.005 
corresponds to 10 km s _1 and v t = 0.2 to 400 km s . 



photons close to line center that are backscattered into the pho- 
tosphere when the resonance zone grows and overlaps with the 
lower boundaryCB Consequently the re-emitted flux in this re- 
gion is higher when calculated via the SEI than when calculated 
via the CMF or MC methods. These discrepancies between the 
CMF and SEI are quite well doc umented and discussed (e.g., 
iHamannll 1 98 U lLamers et al. I ll987h . however we still emphasize 
that one should exercise caution when applying the SEI method 
with high microturbulence on wind resonance lines. Especially 
today, when increased computer-power enables us to compute 
fast solutions using both methods, the CMF is preferable. 

Next we calculated line profiles for structured, ID winds. 
Profiles computed with all three methods agreed for weak and in- 
termediate lines. For strong lines, the agreement between MCS- 
1D and the method from POF, which uses a Sobolev source 
function accounting for multiple-resonance points, was satisfac- 
tory. However, minor discrepancies between Sobolev and non- 
Sobolev treatments occurred for the strong line also when no 
microturbulent velocity was applied (see Fig. IA.2I ). as opposed 
to the smooth case. 

Finally we performed a simple test of our MC-2D code by 
applying it on models in which all lateral slices had the same 
radial structure, i.e., the wind was still spherically symmetric and 
all observers ought to see the same spectrum. We confirmed that 
indeed so was the case, both for smooth and structured models 
(in Fig. lA.2l the latter case is demonstrated). 



10 Remember that neither the SEI nor the CMF, as formulated here, 
include a transition to the photosphere, but treat the lower boundary as 
sharp with a minimum velocity v m j„. 
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Appendix B: The effective escape ratio 

We define the ratio of the velocity gap Av between two clumps 
(see Fig.|6]in the main paper) and the thermal velocity v t as 

Av 

r,= —. (B.l) 

v t 

In the following, we derive an expression for 77, for the wind 
geometry used throughout this paper. If Av tot = Av + \6v\ is the 
velocity difference between two clump centers, we may write 
(omitting the absolute value signs here and in the following) 

Av tot 5v 

Av = Av, ot -Sv= AvtotjS - — 6vp, (B.2) 

Avtot,^ ovp 

where we have normalized the arbitrary velocity intervals to the 
corresponding /3 intervals, suffixes are used to denote parame- 
ters of a smooth velocity law. For notational simplicity we write 

Avtot <5v 
?i = XT*"' ft = T7 . (B.3) 

Assuming radial photons, Av may be approximated by 

Av « Ar totj8 (ft - & —2-), (B.4) 
or Ar tolfi 

with the notations of r following those of v. The volume filling 
factor for the geometry in use is 

/v = tt^- ~ (B.5) 
Vtot r|Ar tot 

with Vci the volume of the clump, V tot the total volume, and 
r\ * r2 the radial points associated with the beginning of 
the clump and the ICM. Using Eq. IB. 51 and Ar tot = vpdt (see 
Sect. 12. 2t , we obtain 

Av * ^vpfittfi - hfv), (B.6) 
or 

and for 77, using the radial Sobolev length of a smooth flow 
U = vj{dvpldr), 

~ t VM hlJV> . (B.7) 



In our models %\ is not given explicitly, but is on the order of 
unity, because we distribute clumps according to the underlying 
smooth /J = 1 velocity law. Thus we approximate 

* (B.8) 



U 

We notice that the porosity length h as defined by 
lOwocki et al.l (l2004t) is h = 1/ f v , where / is the length associ- 
ated with the clump. For the geometry used here this becomes 
h « Sr/f v ~ vpSt. Hence, using ^2 - \ for a smooth velocity 
field, 77 represents the porosity length corrected for the finite size 
of the clump, and divided by the radial Sobolev length. 



